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Abstract 

Several  improvements  have  been  incorporated 
into  the  XAIR  overset-mesh  flow  solver  which  sub¬ 
stantially  improve  the  capability  of  the  code  to  per¬ 
form  time-accurate  solutions.  These  improvements 
include  implementing  an  upwind  flux  formulation 
combined  with  a  quasi-Newton  relaxation  time¬ 
stepping  strategy.  The  code  includes  a  very  stable 
k-e  turbulence  model  which  can  be  used  with  or 
without  wall  function  boundary  conditions.  Results 
are  presented  for  two  moving  body  cases  and  two 
unsteady  flows.  The  new  code  has  demonstrated 
reductions  of  CPU  times  over  traditional  alternating 
direction  implicit  (ADI)  methods  by  factors  of  10  for 
inviscid  flows  and  by  orders  of  magnitude  for  vis¬ 
cous  flows. 

Introduction 

In  order  to  efficiently  compute  the  trajectories  of 
several  bodies  moving  with  mutual  aerodynamic 
influence,  a  flow  solver  must  be  capable  of  taking 
very  large  computational  time  steps  (i.e.,  running 
with  very  large  Courant  numbers).  This  demand  is 
especially  true  when  turbulent  flows  at  flight  Rey¬ 
nolds  numbers  are  computed.  The  time  step  one 
would  like  to  take  in  such  should  be  dictated  by  con¬ 
siderations  of  numerical  accuracy  rather  than 
numerical  stability.  When  an  overlapping  mesh”* 
strategy  is  used,  the  solutions  on  each  mesh  must 
be  tightly  coupled  in  order  to  assure  that  the  effects 
of  the  mutual  aerodynamic  influence  are  properly 
communicated  across  the  interpolated  domain 
boundaries.  With  this  In  mind,  an  effort  sponsored 
by  Dr.  Leonidas  Sakell  of  the  Air  Force  Office  of  Sci¬ 
entific  Research  (AFOSR)  was  initiated  at  the 
Arnold  Engineering  Development  Center  (AEDC)  to 
develop  the  computational  tools  necessary  to  simu- 
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late  the  ripple  launch  of  three  MK-82  stores  carried 
in  tandem  from  an  F-15E  aircraft  at  flight  Reynolds 
numbers  within  a  time  frame  of  a  few  calendar 
weeks.  This  requires  the  solution  of  the  Navier- 
Stokes  equations  with  a  higher-order  turbulence 
model  since  the  aft  two  stores  are  immersed  in  the 
wake  of  the  store  or  stores  in  front  of  them.  Initial 
estimates  of  the  CPU  time  requirements  for  this 
problem  using  the  existing  XAIR^  overset-mesh 
flow  solver  indicated  that  calendar  years  would  be 
required  to  obtain  a  solution  for  the  viscous  ripple 
launch.  Hence,  an  effort  was  made  to  improve  the 
efficiency,  accuracy,  and  numerical  stability  of  the 
XAIR^  overset-mesh  flow  solver.  This  paper  pre¬ 
sents  several  applications  of  the  improved  code  to 
complex  aerodynamic  configurations. 

Code  Improvements 

The  first  improvement  to  the  XAIR  code  con¬ 
sisted  of  replacing  the  central  difference  flux  formu¬ 
lation  with  the  upwind  Harten/Lax/van  Leer/Einfeldt 
(HLLE)  scheme. The  HLLE  scheme  Is  an 
approximate  Godonov  scheme  similar  in  resolution 
to  Roe’s®  flux  difference  scheme.  The  HLLE 
scheme  was  chosen  over  the  Roe  scheme 
because  of  its  lower  operation  count  and  because  it 
allows  for  more  direct  control  of  the  anti-diffusion 
terms  which  are  present  in  both  schemes.  The  lat¬ 
ter  property  is  desirable  when  dealing  with  the  car¬ 
buncle  phenomena  associated  with  both  schemes 
In  hypersonic  shock  capturing-applications.  The 
HLLE  scheme,  like  the  Roe  scheme,  captures 
shocks  crisply  and  has  very  little  numerical  diffu¬ 
sion  in  boundary  layers.  This  makes  these  upwind 
schemes  more  accurate  than  traditional  central  dif¬ 
ference  schemes  with  added  second-  and  fourth- 
order  numerical  smoothing.  Higher-order  spatial 
accuracy  is  achieved  by  using  the  Monotone 
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Upstream  Scalar  Conservation  Law  (MUSCL) 
approach  of  van  Leer.® 

The  second  improvement  involved  replacing 
the  first-order  accurate  in  time  Beam  and 
Warming^  alternating  direction  implicit  (ADI) 
scheme  with  a  second-order  accurate  in  time 
quasi-Newton  approach.  The  Steger- Warming  flux 
linearizations  are  used  to  approximate  the  inviscid 
part  of  the  flux  Jacobians  required  by  the  Newton 
method.  As  shown  by  Whitfield,  et  al.®,  the  use  of 
Steger-Warming  linearizations  provides  the  numer¬ 
ical  stability  required  to  take  large  time  steps.  The 
viscous  part  of  the  Newton  flux  Jacobians  is  pro¬ 
vided  by  linearizing  the  thin-layer  viscous  fluxes. 
Source  term  Jacobian  contributions  for  the  k-e^  and 
Spalart-Allmaras"^®  turbulence  models  are  provided 
by  approximations  of  the  linearized  source  terms. 
The  tridiagonal  system  of  linear  equations  which 
must  be  solved  at  each  iteration  of  the  Newton  pro¬ 
cess  is  treated  using  the  Jacobi  iterative 
method^^’^^  with  underrelaxation.  Second-order 
time  accuracy  is  achieved  by  using  three-point 
backward  differencing. 

All  computational  boundaries  within  the  compu¬ 
tational  domain  are  updated  within  the  Newton 
loop.  This  causes  explicitly  applied  boundary  condi¬ 
tions  to  behave  in  a  pseudo-implicit  manner.  Thus, 
at  the  completion  of  each  time  step,  the  flow  vari¬ 
ables  at  the  boundaries  are  at  the  same  time  level 
as  the  interior  points.  Periodic  boundaries  such  as 
slits  (C-mesh  topology),  overlaps  (0-mesh  topol¬ 
ogy),  and  symmetry  planes  are  treated  implicitly 
through  the  introduction  of  ghost  cells.  This  elimi¬ 
nates  the  “ringing”  at  these  boundaries  often 
encountered  with  explicit  periodic  boundary  condi¬ 
tions  and  allows  the  use  of  higher-order  fluxes  at 
these  boundaries.  Global  coupling  between  the 
overlapping  grids  is  also  handled  within  the  Newton 
loop,  so  that  at  the  end  of  each  time  step  the  entire 
computational  domain  has  advanced  to  the  same 
time  level. 

The  turbulence  models  (k-e^  and  Spalart- 
Allmaras^®)  are  also  Included  within  the  Newton 
loop.  The  Newton  methodology  thus  provides  a 
pseudo-coupling  of  the  turbulence  equations  to  the 
mean  flow  equations.  The  wall  function  boundary 
conditions  of  Nichols®  are  available  for  both  turbu¬ 
lence  models. 


These  improvements  to  the  XAIR  code  have 
substantially  increased  the  allowable  time  step 
which  the  flow  solver  can  take.  The  new  code  has 
demonstrated  the  ability  to  take  time  steps  10-20 
times  larger  than  its  predecessor  for  inviscid  calcu¬ 
lations,  and  is  capable  of  orders  of  magnitude 
improvements  for  viscous  flow.  Local  time  stepping 
is  not  required  for  numerical  stability  with  the  new 
code,  so  all  problems  may  be  treated  as  time-accu¬ 
rate  calculations.  This  has  provided  insight  in  sev¬ 
eral  applications  at  AEDC  which  would  not  have 
been  available  with  a  code  which  requires  local 
time  stepping.  More  details  on  the  new  code  can  be 
found  in  Ref.  13. 

Applications 

Three  unsteady  and  two  moving  body  problems 
are  presented  which  demonstrate  the  capabilities  of 
the  new  flow  solver.  The  unsteady  test  cases  include 
vortex  shedding  from  a  circular  cylinder  at  high  Rey¬ 
nolds  number,  and  a  simple  cavity.  The  moving 
body  problems  are  the  launch  of  a  single  generic 
store  and  the  ripple  launch  of  three  generic  stores. 

The  unsteady  cases  presented  here  utilize  con¬ 
ventional  two-equation  turbulence  models.  The  dif¬ 
ferential  equations  solved  for  these  turbulence 
models  include  the  proper  unsteady  terms,  but  the 
question  still  remains  as  to  how  applicable  these 
turbulence  models  are  to  unsteady  flows.  A  turbu¬ 
lence  model  applicable  to  unsteady  flow  should  be 
capable  of  differentiating  between  unsteady  length 
and  time  scales  and  the  sub-grid  size  turbulence. 
The  turbulence  model  should  only  try  to  simulate 
the  effects  of  the  sub-grid  turbulence.  Conventional 
one-equation  and  two-equation  turbulence  models 
actually  simulate  the  effect  of  turbulent  scales 
larger  than  the  computational  cell  size  and  provide 
no  mechanism  for  filtering  out  the  large-scale  con¬ 
tribution.  Large  Eddy  Simulation  (LES)  turbulence 
models  attempt  to  provide  a  filter  for  the  turbulent 
scales,  but  these  models  have  yet  to  mature  to  the 
point  of  application  for  complex  three-dimensional 
flow  problems.  It  is  expected  that  the  conventional 
turbulence  models  will  overpredict  the  effect  of  tur¬ 
bulence  in  unsteady  applications  and  overdamp  the 
unsteady  nature  of  the  flow.  The  amount  of  this 
overdamping  cannot  be  accurately  determined  at 
this  time,  and  is  probably  dependent  on  the  formu¬ 
lation  of  the  particular  turbulence  model. 
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Vortex  Shedding  From  a  Circular  Cylinder 

The  first  unsteady  flow  test  case  was  the  flow 
past  a  circular  cylinder  with  a  free-stream  Mach 
number  of  0.2  and  a  Reynolds  number  of  8  x  10® 
based  on  the  diameter  of  the  cylinder.  Experimental 
data  for  this  high  Reynolds  number  condition  have 
been  obtained  by  Jones,  et  al.^'^  and  by  Roshko.^® 
This  case  was  chosen  because  the  flow  on  the  cyl¬ 
inder  is  almost  fully  turbulent  and  a  definite  vortex 
shedding  frequency  is  present.  Jones,  et  al.^^, 
identify  three  regimes  for  high  Reynolds  number 
flow  over  a  circular  cylinder  (Re  >  1  x  10®)  in  terms 
of  the  Strouhal  number  (St): 

1.0  X  10®  <  Re  <  3.5  X  10®Wide  band  unsteady 
lift  force  (wide  St  bandwidth) 

3.5  X  1 0® <  Re  <  6.0  X  1 0®  Narrow  band  random 
lift  forces  (narrowing  St  bandwidth) 

6.0  X  10®  <  Re  <  1.0  X  10^  Quasi-periodic  lift 
forces  (small  St  bandwidth) 

Hence  it  is  expected  that  the  lift  force  on  the  cyl¬ 
inder  will  oscillate  at  a  predominant  frequency  for 
this  test  case.  Unfortunately,  the  available  data  also 
indicate  a  sudden  drag  rise  near  this  condition,  and 
large  variations  in  the  drag  coefficient  are  reported. 
The  pressure  taps  used  by  Jones,  et  al.^^  were 
located  only  one  cylinder  diameter  from  the  tunnel 
floor.  Flow  visualization  results  indicate  some  three 
dimensionality  was  indeed  present  in  the  experi¬ 
ment.  The  data  of  Roshko"'®  were  taken  in  a  solid 
wall  tunnel  with  a  model  blockage  (the  ratio  of 
model  frontal  area  to  tunnel  cross-sectional  area)  of 
about  19  percent.  A  10-percent  reduction  was 
made  in  the  drag  coefficient  to  correct  the  data  for 
wall  effects,  but  the  author  admits  that  the  correc¬ 
tions  have  questionable  validity  at  these  conditions 
for  unsteady  force  measurements.  The  data  of 
Jones,  et  al.^'*  were  taken  in  a  slotted  wall  tunnel, 
and  also  had  a  model  blockage  of  about  19  per¬ 
cent.  No  wall  corrections  were  made  for  the  data. 
The  cylinder  used  by  Jones  had  a  smoother  sur¬ 
face  than  did  the  model  of  Roshko.  The  large  model 
blockage  may  contribute  to  the  large  variation  in 
reported  drag  coefficients.  It  would  be  expected 
that  wall  interference  would  cause  the  flow  to  over¬ 
accelerate  at  the  top  and  bottom  of  the  cylinder. 
This  would  cause  a  reduction  in  the  pressure  in  the 


suction  peak  and  probably  delay  flow  separation. 
These  two  effects  would  cause  the  base  pressure 
to  drop,  and  hence  lead  to  a  larger  value  of  drag. 
Thus,  one  would  expect  the  drag  coefficient 
reported  by  Roshko’®  to  be  larger  than  that  of 
Jones,  et  al.’'*,  and  that  both  would  be  larger  than 
an  interference-free  result  as  calculated  here.  As 
will  be  seen,  this  is  the  case.  The  Strouhal  number 
reported  by  both  investigations  is  in  reasonable 
agreement.  Hence,  the  ability  to  predict  the  Strou¬ 
hal  number  will  be  used  as  one  metric  for  evaluat¬ 
ing  the  code. 

Calculations  were  performed  with  a  401  x  201 
grid  with  an  initial  y+  of  50  using  the  k-e  turbulence 
model  with  wall  function  boundary  conditions,  and 
with  a  401  x  21 1  grid  with  an  initial  y+  of  1  using  the 
low  Reynolds  number  k-e  turbulence  model  (inte¬ 
grating  the  turbulence  equations  to  the  wail).  A 
physical  time  step  of  4.3  x  10^  secs  was  used  for 
both  grids,  resulting  in  a  maximum  Courant  number 
of  4770  for  the  wall  function  grid  and  a  Courant 
number  of  2.32  x  10®  for  the  low  Reynolds  number 
grid.  The  calculations  were  begun  from  uniform 
free-stream  conditions.  A  periodic  condition  was 
reached  after  about  1 ,000  iterations,  and  the  solu¬ 
tion  was  advanced  another  2,048  steps  to  provide 
averaged  results. 

Figure  1  shows  the  Mach  number  contours  at  a 
particular  instant  in  time.  The  periodic  disturbance 
is  seen  to  dissipate  as  it  moves  downstream  until  a 
steady  wake  is  finally  achieved.  Unfortunately, 
there  are  no  experimental  data  available  with  which 
to  compare  the  disturbance  dissipation  distance.  It 
is  interesting  to  note  that  Roshko’®  made  hot-wire 
measurements  seven  diameters  behind  the  cylin¬ 
der  and  noted  unsteady  flow  still  present.  The 


Fig.  1 .  Mach  number  contours  for  a  circuiar  cyiin- 
der  at  M  =  0.2  and  Re  =  8  x  1 0®. 
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results  presented  still  show  unsteady  flow  seven 
diameters  behind  the  cylinder.  Roshko  also 
reported  no  unsteady  flow  was  present  at  this  loca¬ 
tion  when  the  Reynolds  number  was  reduced  to 
around  1  x  10®.  In  an  earlier  study,  Delaney  and 
Sorenson^®  discovered  unsteady  flow  at  a  Rey¬ 
nolds  number  of  1  x  1 0®  using  a  hot  wire  located  two 
diameters  behind  the  cylinder.  Hence,  it  seems  that 
the  unsteady  flow  dissipates  rapidly  at  these  high 
Reynolds  numbers  and  the  length  over  which  this 
dissiption  occurs  is  Reynolds  number  dependent. 
Detailed  data  on  the  dissipation  distance  for  the  cyi- 
inder  unsteady  wake  would  be  useful  in  evaluating 
turbulence  models  for  unsteady  flow  applications. 

The  computational  and  experimental  average 
drag  coefficient  and  the  unsteady  lift  coefficient 
Strouhai  number  are  shown  in  Table  1.  The  aver¬ 
age  drag  coefficient  is  underpredicted  by  CFD, 
which  may  be  attributable  to  either  an  overpredic¬ 
tion  of  the  base  pressure  or  a  premature  prediction 
of  the  boundary-iayer  separation  point.  As 
expiained  above,  the  resuits  are  consistent  with  the 
amount  of  waii  interference  present  in  the  experi¬ 
mental  and  computational  studies.  The  Strouhai 
number  predictions  are  in  reasonabie  agreement. 
The  spectrum  of  the  CFD  results  is  shown  in  Fig.  2. 
More  than  90  percent  of  the  energy  is  seen  to  be 
concentrated  in  the  first  spectrai  peak,  indicating  a 
truly  periodic  result. 

Cavity 

A  single  computation  using  the  wall  function 
boundary  condition  in  conjunction  with  the  k-e  tur- 
buience  modei  was  performed  for  a  three-dimen¬ 
sional  cavity.  The  cavity  had  a  length-to-depth  ratio 
of  4.5.  The  flow  computations  were  performed  for  a 
free-stream  Mach  number  of  1 .2  and  a  Reynoids 
number  of  2.0  x  10®  per  foot.  The  characteristic 
time,  tc  (the  time  required  for  free-stream  flow  to 


Table  1.  Computational  and  Experimental  Results 
for  a  Circular  Cylinder  at  M  =  0.2  and 
Re  =8x10® 


Average 

Cd 

Strouhai 

Number 

Wall  Function  CFD 

0.382 

0.293 

Low  Reynolds  Number  k-e  CFD 

0.425 

0.273 

Data  Jones  et  al.''^ 

0.50-0.55 

0.29-0.31 

Data  Roshko  et  al. 

0.68-0.75 

0.27 

traverse  the  cavity),  is  0.00144  secs.  Data  used  for 
comparisons  were  obtained  during  the  Weapons 
Internal  Carriage  Separation  Program  (WICS) 
sponsored  by  the  Wright  Laboratory/Armament 
Directorate.^^  The  basic  configuration  tested, 
shown  in  Fig.  3,  was  a  generic  flat  plate  and  cavity. 
Both  static  and  dynamic  pressure  data  were 
obtained  during  the  test. 

A  coarse  grid  of  approximately  2x10®  points 
was  used  to  model  half  of  the  bay  and  flat  plate.  A 
time  step  of  8  x  10'®  seconds  was  used  for  the  cal¬ 
culations.  The  time  step  is  approximately  50  times 
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Fig.  2.  Spectrum  of  the  normal  force  coefficient 
predictions  for  a  cyiinder  at  M  =  0.2  and 
Re  =  8x  10®. 
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Fig.  3.  Dimensions  of  the  WICS  flat-plate/cavity 
model. 
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larger  than  the  time  step  used  to  calculate  this  flow 
by  Suhs,  et  al.^®  using  the  Ref.  2  flow  solver  with 
the  same  turbulence  model  and  wall  function 
boundary  conditions.  The  calculated  instantaneous 
pressure  on  the  back  wall  of  the  cavity  is  shown  in 
Fig.  4.  The  transients  due  to  the  initial  conditions 
are  seen  to  disappear  after  approximately  20  time 
characteristics,  after  which  the  flow  is  periodic. 
About  0.01  secs  of  computational  data  taken  after 


Fig.  4.  Predicted  pressure  on  the  aft  wall  of  the 
WIC  bay. 


a.  Front  wall 


the  flow  had  become  periodic  was  used  to  provide 
time-averaged  and  spectral  data  for  the  weapons 
bay.  The  averaged  pressure  coefficients  on  the  bay 
centerline  are  shown,  along  with  experimental  data 
in  Fig.  5.  The  computed  and  experimental  sound 
pressure  levels  are  shown  in  Fig.  6.  The  agreement 
is  good  for  both  of  these  quantities  and  consistent 
with  Suhs’  results.  The  spectrum  for  a  point  on  the 
centerline  of  the  bay  back  wall  is  shown  in  Fig.  7, 
along  with  the  first  two  experimental  spectral 
peaks.  It  should  be  noted  that  the  background 
noise  sound  pressure  level  generated  by  the  holes 
of  the  porous  wind  tunnel  walls  is  about  120  dB. 
Again,  the  agreement  is  quite  good.  For  this  case, 
the  total  savings  in  CPU  time  using  the  new  solu¬ 
tion  algorithm  was  a  factor  of  ten  over  that  of  the 
ADI  scheme. 

Generic  Store  Release 

The  first  moving  body  case  chosen  for  presenta¬ 
tion  is  the  generic  wing/pylon/finned  store  configu- 


c.  Back  wall 

Fig.  5.  Time-averaged  pressure  coefficient  for  the  WICS  bay  centerline. 
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a.  Front  wall 


x,ft 

b.  Bottom  wall 


c.  Rear  wall 

Fig.  6.  Sound  pressure  level  for  the  WICS  bay  centerline. 


ration.^®  The  store  physical  properties  and  ejector 
forces  are  given  in  Ref.  1 9.  Euler  calculations  were 
performed  for  a  free-stream  Mach  number  of  0.95 
on  a  1 .5  X  10^  point  grid  system.  The  chimera  mov¬ 
ing  body  methodology  of  Jordan,  et  al.^®  was  used 
to  model  the  store  launch.  A  computational  time 
step  of  0.001 65  secs  was  used  with  the  new  code, 
resulting  in  a  maximum  Courant  number  of  4000  for 
the  grid  system.  This  time  step  is  16  times  larger 
than  that  used  by  previous  investigators. The 
predicted  carriage  loads  are  presented  in  Table  2. 
The  pressure  coefficient  distribution  for  the  store  at 
carriage  is  shown  in  Fig.  8.  Results  of  the  computed 


Frequency,  Hz 


Fig.  7.  Spectrum  of  the  sound  pressure  level  on 
the  back  wall  of  the  WIGS  bay. 


and  experimental  trajectories  are  jg|j|Q  2.  Predicted  and  Measured  Carriage  Loads  for  the  Single 
shown  in  Figs.  9-1 2.  In  all  cases,  the  Generic  Store 


agreement  is  excellent.  The  results 
are  consistent  with  those  of  Lijewski^® 
and  Thoms  and  Jordan  For  this 
case,  the  total  savings  in  CPU  time 
was  a  factor  of  six  over  that  reported 
by  Thoms  and  Jordan.^^ 
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TER  Ripple  Launch 

The  final  case  presented  is  a  multipie-body 
reiease  from  a  tripie  ejector  rack  (TER)  configura¬ 
tion.  The  case  included  the  wing/pylon  configura¬ 
tion  from  the  previous  test  case,  along  with  three 
generic  stores.  The  store  physical  properties  and 


Fig.  8.  Pressure  coefficient  contours  on  the 
generic  store  at  carriage. 


Fig.  9.  Comparison  of  computed  and  measured 
force  coefficients  for  the  single  generic 
store. 


ejector  forces  are  given  in  Ref.  21.  The  bottom 
store  was  reieased  first  at  t  =  0.000  sec,  followed  by 
the  outboard  store  at  t  =  0.040  sec,  and  finally  the 
inboard  store  at  t  =  0.080  sec.  The  grid  system 
included  3.2  x  10®  points.  The  flight  condition 
selected  for  this  demonstration  was  a  free-stream 
Mach  number  of  0.95  and  an  altitude  of  20,000  ft. 
The  k-e  turbulence  model  was  used,  along  with  the 
wall  function  boundary  conditions.  The  computa¬ 
tions  were  performed  using  a  time  step  of  0.0033 
secs,  which  is  a  factor  of  32  times  larger  than  the 
time  step  used  in  the  Euler  calculations  of  Thoms 
and  Jordan The  maximum  Courant  number  was 
1  X  1 0'*.  Pressure  coefficient  distributions  with  the 
stores  at  carriage  are  shown  in  Fig.  13.  The  pre¬ 
dicted  carriage  loads  are  shown  in  Table  3.  The 
present  viscous  loads  are  generally  smaller  than 
the  Euler  results  of  Thoms  and  Jordan.^*  This  is 
consistent  with  the  fonward  movement  of  the  shock 
in  the  tail  region  of  the  stores  due  to  thickening  and 
separation  of  the  boundary  layer.  The  computa¬ 
tional  predictions  for  the  inboard,  outboard,  and 
bottom  store  trajectories  and  Euler  angles  are 


Fig.  11.  Comparison  of  computed  and  measured 
store  trajectory  for  the  single  generic 
store. 


c  1.5 


0.5 
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1  -0-5 
£ 

S  -1.5 


-2.5 


I  I  I  I 

Roll  Euler  CFD 

Roil  Exp - 

Pitch  Euler  CFD  >« 
Pitch  Exp.  ■■■■■■ 
Yaw  Euler  CFD 
Yaw  Exp. 

. 


0.05  0.10  0.15  0.20  0.25  0.30  0.35 
Time,  sec 


Fig.  10.  Comparison  of  computed  and  measured 
moment  coefficients  for  the  single  generic 
store. 


Fig.  12,  Comparison  of  computed  and  measured 
store  orientation  for  the  single  generic 
store. 


7 

American  Institute  of  Aeronautics  and  Astronautics 


shown  in  Figs.  14-19.  The  predicted  roll  angle  for 
the  bottom  store  at  the  end  of  the  trajectory  is  about 
three  times  larger  than  that  predicted  by  Thoms  and 
Jordan.^^  The  predicted  roll  angle  of  the  inboard 
and  outboard  stores  is  opposite  of  that  predicted  by 
Thoms  and  Jordan.^^  The  pitch  angle  results  are 
similar  between  the  two  predictions  as  expected 
since  they  are  dominated  by  the  ejectors.  It  should 
be  noted  that  the  present  Navier-Stokes  calcula¬ 
tions  required  one-third  of  the  CPU  time  required 
for  the  Euler  predictions  of  Thoms  and  Jordan.^^ 


Fig.  13.  Pressure  coefficient  contours  on  three 
generic  stores  at  in  TER  configuration. 


Conclusions 

Several  improvements  have  been  incorporated 
into  the  XAIR  overset  mesh  flow  solver  which  sub¬ 
stantially  improve  the  capability  of  the  code  to  per¬ 
form  time-accurate  solutions.  These  improvements 
include  implementing  an  upwind  flux  formulation 
combined  with  a  quasi-Newton  relaxation  time¬ 
stepping  strategy.  In  addition,  use  is  made  of  the 
wall  function  boundary  condition  for  both  the 
Spalart-Allmaras  and  k-e  turbulence  models  to 
reduce  the  number  of  grid  points  in  the  boundary 
layer  and  to  reduce  the  maximum  Courant  number 
occurring  in  a  computational  mesh.  The  new  code 
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Fig.  15.  TER  inboard  store  orientation. 
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Fig.  14.  TER  inboard  store  trajectory. 


Fig.  16.  TER  outboard  store  trajectoiy. 


Table  3.  Predicted  Carriage  Loads  for  the  Stores  in  a  TER  Configuration 


CA 

CY 

CN 

Cl 

Cm 

Cn 

1.134 

^.266 

0.286 

0.001 

0.535 

Bottom  Store  Present 

0.974 

-0.243 

0.137 

-0.023 

0.537 

Inner  Store  Thoms^^ 

1.382 

-0.704 

1.093 

-0.060 

-2.999 

1.094 

Inner  Store  Present 

0.945 

-0.626 

0.709 

0.007 

-2.047 

Outer  Store  Thoms^^ 

1.180 

-0.432 

0.118 

-0.010 

-1.001 

0.911 

Outer  Present 

1.007 

-0.313 

-0.052 

-0.004 

-0.723 

0.789 
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has  been  tested  on  a  variety  of  steady  and 
unsteady  flows.  Four  of  these  test  cases  are  pre¬ 
sented.  The  improvements  have  allowed  computa¬ 
tions  to  be  performed  with  time  steps  ten  times 
larger  than  previous  ADI  methods  for  Euler  flows, 
and  has  increased  the  maximum  time  step  by 
orders  of  magnitude  for  viscous  flows.  Further  vali¬ 
dation  efforts  are  underway  and  will  be  reported  at 
a  later  date. 
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